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Abstract 

We present a flexible Monte Carlo implementation of the perturbative framework of High 
Energy Jets, describing multi-jet events at hadron colliders. The description includes a 
resummation which ensures leading logarithmic accuracy for large invariant mass between 
jets, and is matched to tree-level accuracy for multiplicities up to 4 jets. The resummation 
includes all-order hard corrections, which become important for increasing centre-of-mass 
energy of the hadronic collision. 

We discuss observables relevant for confronting the perturbative framework with 7 TeV 
data from the LHC, and the impact of the perturbative corrections on several dijet and trijet 
observables which are relevant in the search for new physics. 
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1 Introduction 

The cross-section at the LHC for particles charged under QCD will generally be larger than 
that for colourless particles, and so many of the discovery channels used in the search for 
new physics involve the detection of hard, hadronic jets. The large mass hierarchy between 
any (often heavy, in order to avoid existing exclusion limits) new particle produced and those 
of the decay products often implies that many jets will be produced in the decay of a new 
state. The finger prints of any such new physics will, however, have to be found amongst a 
large contribution to the same signature from multi-jet processes within the Standard Model. 
Therefore, a detailed understanding of the Standard Model processes will assist in the search 
for new physics. Examples of Standard Model processes acting as background to many searches 
for new physics are e.g. W, Z-|-jets (especially with 3,4 jets or more). 

However, even the nature of some Standard Model processes is best studied in events with 
multiple jets. For example, the CP-structure of the induced Higgs boson couplings to gluons 
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through a top-loop could be measured by a study of the azimuthal angle between the two jets 
in events with a Higgs boson in association with dijets [HE]- 

In both examples, hard radiative corrections will be sizeable at the LHC, by which we mean 
that the exclusive (n + l)-jet rate is a significant component of the inclusive n-jet rate. And 
more so in many of the regions of interest for searches of new physics. Therefore, a tree-level 
description of the inclusive re-jet process will be unsatisfactory for the involved analyses beyond 
just a measurement of the cross section. 

The reason for the increased importance in many situations of hard, perturbative corrections 
at the LHC over the situation at previously, lower energy colliders is very simple. Two effects 
act to suppress hard corrections: the increasing powers of the perturbative coupling, and the 
increase in the light-cone momentum fraction of the partons extracted from the proton beyond 
that necessary for the final state without the additional hard jet. The suppression from this last 
kinematic effect is caused by the decrease in the parton density functions (pdf) as the light-cone 
momentum fraction x is increased. However, for processes with at least two particles in the final 
state, there is a fine trade-off between the suppression from the pdf and the increasing phase 
space for additional emission in-between the most forward/backward hard jet (even when this 
additional emission is hard in transverse momentum), as the rapidity span between the two most 
forward/backward jets is increased. At previous, lower-energy colliders, this balance was tipped 
more towards a suppression than will be the case at the LHC. 

At previous colliders, the "significant" rapidity separation of the two objects, which is neces- 
sary for the opening of phase space for additional radiation, would already bring the light-cone 
momentum fractions into the region of extremely fast falling pdfs as x — )• 1, thus effectively ve- 
toing additional emissions. However, the situation is different for the LHC processes discussed 
above, since in the case of e.g. M^-boson production with at least 3 jets, two jets will naturally be 
produced with a size-able separation in rapidity [3]. In this case, there is only a small suppres- 
sion for additional (especially central) radiation, even when the additional jets have a sizeable 
transverse momentum. This holds true also for other processes, provided the hard scattering 
amplitude has a mechanism for effectively radiating into the rapidity span. This is the case 
when colour is exchanged between the particles either side of the span, whereas a colour-singlet 
exchange leads to less radiation in the span [ll[5]. 

While a fixed order (e.g. LO or NLO) calculation may be adequate for the description of 
sufficiently inclusive quantities like the total inclusive cross section, the question is to what ex- 
tent a given theoretical description allows for the radiation into the phase space which becomes 
available with the increase in partonic centre of mass energy — and how important the descrip- 
tion of this radiation is for a given observable. The current paper discusses these problems, and 
presents results obtained in a recently proposed all-order perturbative framework. 

It is clear that a NLO calculation allows for just one, also hard, additional emission above 
the minimum number of jets required in the analysis. The all-order description of a parton 
shower, on the other hand, captures the soft and collinear emissions, but will underestimate the 
amount of hard radiation. This deficiency can be repaired order-by-order through a CKKW-L- 
style matching [6HH], or to full NLO accuracy [HllIO] for low multiplicities. In both cases, the 
deficiency of the parton shower in describing hard radiation is repaired by the use of full tree- 
level matrix elements. The maximum multiplicity applied in the tree-level matching is limited 
by the time for evaluation of the tree-level matrix elements. Since in a CKKW-L-style analysis 
the matching scale should be chosen somewhat smaller than the transverse scale required in the 
definition of jets to avoid matching artefacts, the matching procedure will run out of available 
matrix elements at a lower multiplicity than the maximum for which the LO process has been 
calculated. This can be viewed simply as a result of the attempt within the matching procedure 



3 



to describe not just the total rate, but also the final state configuration. 

The framework of High Energy Jets (HEJ) |liyi2j provides an all-order description of pro- 
cesses with more than two hard jets, based on an approximation which captures the hard, 
wide-angle emissions missed in a shower-approach based on soft and collinear splitting func- 
tions. HEJ does not try to redo the job of the shower, but focuses specifically on the part 
not done by a parton shower. Work is in progress to combine the description of HEJ with a 
parton shower [13J; the most important component of the matching between HEJ and a parton 
shower is the avoidance of double counting of soft radiation, which is treated to all orders in 
both descriptions. 

The formalism of HEJ is inspired by that underlying |14H16j the BFKL equation |17j , and as 
such, an approximation for both real and virtual corrections is obtained to all orders, obviously 
with all IR divergences cancelling between the two contributions. Differently to the BFKL 
approach, however, HEJ applies an approximation only to the partonic scattering amplitudes, 
and not the phase space integration, which is performed for each explicit multiplicity. In this 
respect, HEJ resembles a parton shower formulation of an all-order summation. Furthermore, 
by applying the approximations at the level of the scattering amplitude Ad (and not |A4p), it is 
possible to supplement |lipi2j the approximations with the requirement of e.g. gauge invariance, 
and thereby obtain a formalism, which reproduces more accurately the fixed order perturbative 
results when checked order by order, while simultaneously being sufficiently simple that all-order 
results can be explicitly obtained. 

In the current paper, we develop further the formalism of High Energy Jets by matching to 
fixed order results and include some sub-leading corrections. Furthermore, we demonstrate the 
application of HEJ to the production of at least two and at least three jets. 

The structure of the paper is as follows: in Section [2] we briefly review the formalism within 
High Energy Jets, which allows approximate all-order results to be obtained |lHll2j. In Sec- 
tion [3 we describe the matching of these amplitudes to full, high-multiplicity tree-level results. 
In Section|3]we include some sub-leading corrections, which stabilise the dependence on the scale 
choice [18iil9j. In Section [5] we present results for dijet production obtained with the full formal- 
ism of High Energy Jets, and discuss observables and distributions for which the higher-order 
corrections are particularly important in order to obtain a perturbatively stable description. 
These can lead to a direct experimental test of the importance of the correct perturbative de- 
scription. 

The all-order results presented in this paper are obtained using the implementation of the 
formalism of High Energy Jets in a fully flexible parton-level Monte Carlo generator, which can 
be downloaded at http://cern.ch/hej. 

2 All Orders with High Energy Jets 

The all-order perturbative framework of High Energy Jets (HEJ) initiated in Ref. |lipi2ll20p21j is 
addressing some of the short-comings in the description of multiple hard, perturbative corrections 
in both the (low) fixed-order and in the parton shower formulation. The perturbative description 
obtained with HEJ reproduces the correct, all-order, full QCD limit for both real and virtual 
corrections to the hard perturbative matrix element for the hard, wide-angle emissions which 
underpin the perturbative description of the formation of additional jets. The central parts of 
the formalism were presented in Ref. jlipi2j and discussed further in Ref. [22p23j . In this section, 
we will give just a brief overview of the formalism on which the approximations are based; the 
next section will then discuss how to incorporate matching corrections to full, high multiplicity 
tree-level accuracy. 
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2.1 Dominance of the t-channel poles, and current-current scatting 

In the standard parton shower formahsm, the physical picture is one of successive branchings off 
s-channel propagators, governed by the DGLAP sphtting functions [24H27j . Such a framework 
can sum to leading logarithmic accuracy, and to all orders the behaviour dictated by the soft 
and coUinear s-channel singularities arising in the perturbative corrections to a given scattering 
amplitude. It describes correctly emissions with small invariant mass to the hard scattering 
amplitude. 

The limit of pure A^-jet amplitudes for large invariant mass between each parton of similar 
transverse momentum is described by the FKL-amplitudes |14pi5]. which are at the foundation 
of the BFKL framework [17j . The physical picture arising from the FKL amplitudes is one of 
effective vertices connected by t-channel propagators. The reduction of the formalism to the 
two-dimensional BFKL integral equation relies on many kinematical approximations, which are 
extended to all of phase space. Using an explicit (or so-called iterative) solution to the BFKL 
equation |28H30] , it is however straightforward to show that despite the logarithmic accuracy (in 
s/t), the perturbative expansion of the (B)FKL solution does not give a satisfactory description 
of the results obtained order by order with the true perturbative series from QCD [21]. 

High Energy Jets pTflSj inherits the idea of effective vertices connected by t-channel currents 
in order to reproduce the correct limit of N-iet amplitudes, but goes beyond controlling just 
the logarithmic accuracy like the FKL formalism. The kinematic building blocks of the FKL 
formalism depend on transverse momenta only, as a result of the kinematic limits applied in order 
to separate the amplitude into effective vertices separated by t-channel exchanges [31] . In the 
following, we will discuss how to obtain a better approximation for the t-channel singularities. 

The 2—7-2 scattering qQ — t- qQ obviously proceeds through just a t-channel exchange of the 
gluon current generated by a quark. A good formalism for the description the t-channel poles 
should get at least this very simple process exact. The colour and helicity averaged and summed 
square of this simple scattering amplitude is given by 

i-r-rtree i2 4 ^ '^'^ ~l~ /-i \ 

\MqQ^,Q\ =9' (1) 

Despite its simplicity, this amplitude can already be used to illustrate the problem of the ap- 
proximations made in the standard BFKL procedure. The limit of Multi- Reg ge- Kinematics 
is defined for the scattering process pa,Pb PIt ■ ■ ,Pn in terms of transverse momenta and 
rapidities y = ln( j as the following conditions 



or alternatively 



Vi E {2, . . . ,n - 1} : yi-i '>yi'> yi+i, ^i,j : \pi±\ ^ \pjj_\, 



Vi, j : Iml ~ \Pj±\, Sij 00, (3) 



where Sij = 2 Pi.pj and s = 2 pa-Pb- For the 2 — )• 2 process, the MRK limit of the Mandelstam 
variables is given by f — )■ —k'j_, s ~ —u — t- oo. The effective approximation applied in the BFKL 

formalism (both at LL and NLL) for the 2 — > 2 process is \M^Q^gQ^'^^\'^ = 5^ 1 •^^/(^i)^- 
However, for much of the kinematics relevant at the LHC, t and —k'j_ differ by at least an order 
of magnitude, and s and u differ significantly, leading to a gross overestimation of the cross 
section, if the BFKL approximation is applied. 

In Eq. (pP), the s^-component arises from scattering of quarks of the same helicities (e.g. q~Q~ - 
q~Q~), whereas the u^-component arises from the scatting of unlike helicities (e.g. q~Q^ — s- 
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q^Q^)- Since this difference is important in obtaining sufficient accuracy, HE J is based on 
the calculation of scattering processes at the amplitude level (as opposed to the square of the 
amplitude), and the sum over helicities is performed explicitly. For the gQ-process then, the 
obvious choice of formalism is that of current-current scattering. 

In the spinor notation for the quark currents (see Ref. [llj for details), = u^^^u^ is 
written as (l|^|a), and then the (colour and coupling stripped) matrix element for the process 
%„QpH 1v^Qv2 reads 



{l\^,\ar—{2\v\h). 



(4) 



While it is possible to shorten this expression by use of the Fierz identity, we choose to keep 
the formulation in terms of currents, as this will prove useful for the generalisation to other 
processes, including W, H, Z+jets. 

Let us denote the spinor string (for helicities ha,hi,hb,h2 of the quarks) appearing in the 
amplitude as 



qhahi,^hlh2 

'^qQ-fqQ 



(1 hi\fl\a ha) g^" (2 h2Hbhb) 



(5) 



This complex number can be calculated using any explicit representation for the spinors (see 
e.g. Ref. [TT|fT2] ). and we will denote the sum over helicities of the absolute square of this number 
by 



\S, 



qQ->qQ\ 



ha,ha,hi,,h2 



qhahij^hlh-z 



(6) 



Of course in this case non-zero contributions arise only when ha = hi and hh = h2- 

The colour and helicity summed and averaged matrix element for the scattering process 
qQ ~^ qQ is then 



l-^qQ-^qQl 



1 



\\SqQ^qQ\ 



*2 



(7) 



with ti = {pa — Pi)^ and t2 = {—pb +^2)^ (obviously ti = t2 in this case of a 2 — )• 2-process), 
which equals Eq. ([1]). 

The point of this tour de force through the simple formalism of gQ-scattering is that using 
this formalism, the amplitudes for g(7-scattering can be recast in a very similar form. In fact, a 
careful analysis [12] of the helicity structure in qg — t- gg-scattering reveals that all the amplitudes 
where the helicity of the gluon is unchangecll] factorise again into two currents contracted over 
a t-channel pole. For example, the fully colour-dressed scattering amplitude for the process 
Q~{Pa) + g^{Pb) q'ivi) + 9^{P2) equals [12] 



M 



■2 P2± { f2 fb iPb ^.b +2 IP2 
\P2l\ \ V P2 



teattexl— 1 ^ ^^l^'^)' (8) 

Pb I * 



^AU helicity-flip amplitudes are systematically suppressed by a factor of s. 
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with = Px + i Py, = E — pz- We have taken the negative z-direction to be that of the 
incoming gluon, without loss of generahty. We immediately recognise the kinematic structure 
(in terms of currents) of g~Q~'"-scattering, multiplied by a momentum-dependent colour factor. 
The colour summed and averaged scattering matrix element is 



where z = P2 /p^ (and again ti = {pa — pi)'^ = {pb — ^2)^ = ^2)- This has a striking similarity 
to the amplitude for qQ-scattering (see Eq. ([7])). In fact, it differs only by the slightly more 
complicated colour factor in square brackets, which replace the Cp in the case of quark scattering. 
In the MRK limit — )• p'^ this tends to Ca, and the qg scattering matrix element is equal to 
the one for qQ rescaled by Ca/Cf, just as expected [32|. Eq. Q is, however, the exact result, 
and the square bracket is strictly larger than Ca, and uniformly decreasing for increasing z. 
Small z here reflects a large change in light-cone momentum for the gluon, and unsurprisingly 
a strong acceleration is reflected in an effectively stronger interaction (though this is unrelated 
to the higher order perturbative effect of the running of the coupling) . 

The perhaps most interesting result of using the formalism of currents directly is the obvious 
display that this process has just a t-channel pole (i.e. no poles in the s or n-channel), exactly 
like the seemingly simpler qQ-scattering process. The pure f-channel structure, and the same 
colour factors, hold true for all the helicity assignments, where the gluon does not flip helicity 
during the scattering. Not only are the helicity-flipping amplitudes systematically suppressed 
by the centre-of-mass energy, half of these amplitudes even still have only a t-channel pole, s 
and u channel poles appear only in the helicity scatterings where the helicity of the incoming q 
and g is identical, and the helicity of the gluon is flipped under the scattering. 

In the case of pure gluon scattering, it could seem a little arbitrary to discuss the s,t, 
and u-channels. However, in the cases of scattering of two gluons of opposite helicities, like 
g^ g^ —7- g^g^ it turns out again that the scattering amplitude has just a t-channel pole, 
and is again just the contraction of two currents with special colour factors, which depend only 
on the acceleration of each gluon during the scattering. 

The brief summary presented here of the studies in Ref. |l H ll2 [ [22 t l23j illustrates how the 
t-channel exchange is completely well-defined (for 2—7-2 processes) not just for gQ-scattering, 
but also for scattering processes involving gluons. This is displayed in a formalism based on 
helicity amplitudes and currents, without resorting to kinematic approximations or limits. 

2.2 All Orders Real Corrections 

The previous section demonstrated that the t-channel pole of the full scattering amplitudes is 
much more important for the accurate description of the scattering processes than the zoo of 
Feynman diagrams would suggest. We demonstrated how this t-channel pole can be described 
exactly for many processes by a formalism based on the scattering of specific helicity currents. 
For example, the colour summed and averaged amplitudes for both qQ — )• qQ and the dominant 
ones of qg — )• qg scattering (all helicity non-flipping, like q~ g'^ — )• q^g^) are described exactly 
by a formalism of pure quark current scatterings, with colour factors depending on the flavour 
(quark or gluon) of the scattered partons. 




(9) 
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In the current section we will describe briefly the approximations to the real, radiative 
corrections of the 2 — )• 2 process in High Energy Jets. The soft and collinear regions are already 
well understood by the description in a parton shower. HE J focuses instead on the hard, radiative 
corrections. The aim is to build a framework which is sufficiently accurate for a "first guess" for 
the impact of the radiative corrections (i.e. to all order with a certain logarithmic accuracy), but 
which then is also sufficiently flexible to include matching to the full fixed-order result, where 
this is accessible. The control of the cross section to leading logarithmic accuracy in log(s/t) 
requires control of the hard scattering matrix element to leading power in s/t, as s/t — t- oo. 
As discussed in the previous section, and in more depth in e.g. Ref. |1 HI121[2T] . the control of 
the leading power alone is achieved already in the formalism of Fadin-Kuraev-Lipatov [14fll6] . 
but this is insufficient to ensure a good description of the scattering amplitudes in the energy 
regime of the LHC. The formalism described here combines the right limit at s/t — )• c« (or 
more generally the MRK limit) with complete gauge-invariance. It is already clear that the 
description of the 2 — )• 2-processes discussed in the previous section is gauge invariant, since it 
describes the t-channel pole of the full scattering amplitude exactly. In the current section we 
will build a gauge-invariant approximation to 2 — n-processes. 

2.2.1 Dominant n-jet configurations 

First, we will discuss briefly which processes dominate the 2 — t- n partonic scattering in the 
MRK limit. For any 2 ^ n scattering process, the final state particles can obviously be ordered 
according to rapidity. Apart from exceptional phase space points (of zero measure), no two 
particles will have the same rapidity. 

At the currently implemented accuracy, the HEJ amplitudes will describe the leading con- 
tribution (in the invariant mass between two neighbouring partons) to the n-jet production 
process. For a given n-jet kinematic configuration, many of the possible partonic channels will 
be systematically suppressed. These channels will not be summed to all orders, but will be 
included "only" through matching corrections. Consider now the rapidity ordered final state 
jets. The leading contributions to the n-jet configurations are those where the fiavour of the 
most forward jet equals that of the incoming parton of positive light-cone momentum, and the 
flavour of the most backward jet is identical to that of the incoming parton of negative light-cone 
momentum. The leading contribution to jet production between the jets extremal in rapidity is 
given by pure gluon emissions. Such processes can proceed through a gluon exchange between 
all rapidity-ordered particles. Changing the flavours of two jets, such that a single gluon propa- 
gator between the two jets is replaced by a i-channel quark propagator, automatically leads to 
a suppression of 1/sij for Sij — )• oo, where Sij is the invariant mass between the two jets. We 
choose to call the leading configurations Fii'L-configurations, since they are the same as those 
which are considered in the amplitudes by Fadin-Kuraev-Lipatov |14H161[33Vt35j . 

2.2.2 Amplitudes and Effective Vertices 

The all-order approximations of the n-parton FKL-configurations are constructed similarly to 
the 2—7-2 scattering amplitudes considered in Section 12. H as effective vertices connected by t- 
channel propagators. In the case of the partons of largest or smallest rapidity, these are directly 
the effective currents discussed in the previous section. The emission of additional gluons is 
performed by gauge-invarianllf], effective vertices. These were derived in Ref. [llj, and take 

^by which we of course mean fully gauge invariant, not just up to sub-asymptotic terms as it is often meant 
in the BFKL literature. 
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into account the leading contribution from emissions off both the t-channel exchange and the 
two incoming and the most forward/backward outgoing partons. The effective vertex for the 
emission of a gluon of momentum Pg = qi — q2, V^iQi, Qi+i), is given by [11] 



VP{qi,qi- 



i+1 



r 



p 

Pb 



_^ Pi+l ■ PB _^ Pi+1 



Pr. 



Pi+l -pA PA- PB 



J2 

ii+1 



^ Pi+1 ■ PA ^ 



PA -Pn 

Pi+1 ■ Pi 



+ PA^Pl 



PB^Pn 



(10) 



2 \Pi+i -pB PB- PA PB ■ Pi 

This form of the effective vertex is fully gauge invariant; the Ward Identity, pg ■ V = Q can 
easily be checked. This allows for a meaningful approximation to the scattering amplitude to 
be constructed. 

Another approximation of HEJ is then the systematic omission of interference effects between 
identical particles, since such effects are suppressed by the invariant mass between the particles. 
Essentially, each emission is treated as a distinguishable particle, just like in a parton shower. 
The resulting tree-level approximation for a 2 — )• n scattering is illustrated in Fig. [TJ Virtual 
corrections modify the t-channel propagators and are discussed together with regularisation in 
the next section. The tree-level HEJ-approximation for the square of the amplitude describing 
a (/Q-scattering process with n jets in the final state is then given by [11] 



M 



qQ^qg...gQ 



4(iV2 



1 



Cf 



1 



(11) 



n 



-g'CA 
tiU+i 



V^{qi,qi+i)V^{qi,qi+i] 



where ||5qQ_s.gQ|p indicates the square of pure current-current scattering of Sec. 12. li In the 
case of scattering of gluons, the terms in this sum are weighted with helicity-dependent colour 
factors [12j, one of which is listed in Eq. Q. All the building blocks for constructing the High 
Energy Jeis-scattering amplitudes are listed in Appendix lAl 

2.3 All Orders Virtual Corrections 

The virtual corrections are approximated with the Lipatov ansatz for the f-channel gluon prop- 
agators (see Ref. [TTJ for more details). The Lipatov Ansatz states that order by order, the 



Pa — Current — Pi 

qi I 



q2 i 



P2 
P3 

Pn-1 



qn-1 

Pb — Current I — Pn 



Decreasing 
Rapidity 



Figure 1: The analytic structure of a scattering amplitude in High Energy Jets. 
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leading logarithmically virtual corrections to the full n-parton scattering amplitude in the MRK 
limit can be obtained by the following replacement in the scattering amplitudes: 

^ ^ exi)[a{qi){yi-i-yi)] (12) 

with 

«('?0 = VQ^|y=#'(qV^T- (13) 

This ansatz for the exponentiation of the virtual corrections in the appropriate limit of the n- 
parton scattering amplitude has been proved to even the sub-leading level [3 HI36H38] . In Sec. H] 
we will discuss parts of the next-to-leading logarithmic corrections, which can be included as 
corrections of the type /3o log((5^//i^). 

2.4 Generation and Regularisation of the Cross Section 

We will now discuss the construction of the all-order, regularised dijet cross section. The neces- 
sary details for constructing a generator were already given in Ref. [II112I], but the discussion 
here is more detailed. We begin by defining the matrix element squared built from the t-channel 
factorised picture (eq. (jlip ) combined with the virtual corrections discussed in the previous 
section: 

2 



■^"^^e /i/2^/i9-9/2 



|2 



1 ||o 

4(A0-1) ll^/i/2->/i/2ii 

"-2 /_„2^^ \ (14) 

. Jl M— i y^(g,, gi+i)V;.((?„ g.+i) 

n-l 

■ Y[ exp [2a{qj){yj-i - yj)] , 
i=i 

where /i, /2 indicate the flavour (quarks or gluon), Sf^f^^f^f^ is the sum of contracted currents, 
and Kf^ is Cp if /i = g and Ca if fi = g- These pieces are all given explicitly in Appendix Rl 

The dijet inclusive cross section is simply constructed as the explicit phase space integral over 
the explicit sum of real, radiative corrections, including the leading, all-order virtual corrections. 
We illustrate the procedure with gQ-scattering, but the generalisation to incoming gluons is 
straightforward using the gluon currents and factors detailed in Appendix Rl 

tv 2 

n=2 i=l \J \ ''^) / ^ ^^^-^ 

X {27T)U'(^pPk?j 02,(fe}). 

Here, {pi±,yi) denote the transverse momentum and rapidity of the i'th final state parton. The 
parton momenta fractions are given by 

Xa = Y] exp(-yi), Xb = Y] exp(yj), (16) 
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with ^/s the total hadronic centre-of-mass energy. In Eq. (jlSp . VI denotes the total partonic 
centre-of-mass energy, s = XaXbS, and fA,q{xa,Qa),fB,QixB,QB) denote the relevant parton 
density functions for parton A, B respectively at the resolution scales Qa,Qb- We will discuss 
the choices of scales further in Sec. |H The function 02j{{pi}) takes as arguments all the final 
state partons, and returns 1 if there are at least two jets, according to the chosen jet-definition. It 
is otherwise zero. In the current study, we choose to apply the anti-kt algorithm as implemented 
in Fast Jet [39], with a i?-parameter of 0.6; however, obviously any jet-definition can be applied 
on the partonic ensemble. 

The integration over transverse momentum runs from to infinity. We choose to generate 
only the rapidity ordered phase space (i.e. yi-i < yi < yi+i) using the approach of Ref. [IQ], since 

the i?i?J-amplitudes jA^e'^'l take as argument the rapidity ordered set {pi}- The phase space 
integration of standard fixed-order amplitudes can be done in a similar way (and indeed is done in 
the matching-procedure of SecE]), where then an additional Monte Carlo sampling is performed 
over the identification between the particle leg and the rapidity ordered set of momenta. The 
phase space generation method of Ref. [ID] is very efficient for processes dominated by t-channel 
poles. 

The matrix elements |A^e'^| are divergent for any pi± — )■ 0. We will first discuss how for all 
but the extremal partons, this divergence cancels with the pole in e from the virtual corrections 
implemented according to the Lipatov Ansatz for the resummed i-channel propagators (we will 
then return to the case of the extremal partons below). 

Consider the limit where the transverse momentum of the iih. emitted gluon is vanishing. In 
this limit. 



■'*'l£ Pa Pb^Pl ■■■ Pi~l Pi Pi+1 ■■■ Pn 



Ca 



■'^'■e Pa Pb -> Pi ■■■ Pi-1 Pi+1 ■■■ Pn 



(17) 



where the matrix element on the RHS has n — 1 final state particles, and p? is the sum of the 
squares of the transverse components of pi in the Euclidean metric. By integrating over the soft 
region p? < of phase space in D = A + 2e dimensions we find 



^d2+2^pdyi [Ag'CA 



(27r)2+2£ 47r V P^ 



;Ayi„i,i+i— — — -(A //i ) . 



(18) 



(2vr)2+2e47r r(l + e) e 



The square of the matrix element on the left hand side of Eq. (jl7p contains the exponential 
exp(2d;(gj)Ayj_i^j4.i). By expanding the exponential to first order in and in e, the resulting 
pole in £ does indeed cancel that of Eq. (jlSp . and the combined effect of one soft real emission 
and the first term in the expansion of the Reggeised propagator is a factor 

A„:_„,.^l„(|) (19) 

multiplying the (n — l)-particle matrix element. It is clear that the nested rapidity integrals of 
additional soft radiation in the t-channel factorised multi-parton amplitudes will build up the 
exponential needed to cancel the poles from the virtual corrections to all orders in q<j. The 
divergence arising from a given real emission is therefore cancelled by that arising from the 
virtual corrections in the Reggeised t-channel propagator of the matrix element without the 
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relevant real emission. Therefore, if indeed Eq. (jl7p had been an equality for pf < A^, then the 
regularised HE J matrix element squared would be: 



n-1 

• JJ exp [a;°(gj, A)(yj_i - y^)\ , 
i=i 

^ {(IjA) = log 



(20) 



vr 



which should only be evaluated for p? > A^, and a simple phase-space slicing would then have 
been sufficient to organise the cancellation of divergences. However, while Eq. (|17|) does describe 
the divergence in the soft limit, it is not an exact identity. We can account for the finite difference 
by including an integration over 

-^V^{qu qi+i)V^{qi, Qi+i) - \ (21) 

for p? < A^. Numerically, it turns out to be sufficient to account for the difference and include 
this integral for values of |pi| above roughly 0.2 GeV. The regulated matrix elements for HEJ 
are then given by 



|-^HEj(te})| = 4 (Ag _ 1) 11^/1/2^/1/211^ 



i i . A-,, 



n—2 / / 1 A \\ 

•11(5'^^ (T^^^fe,'?^+i)^M(9^,g^+i)- — ^(P?<A'))) (22) 

n-1 

• exp [u}^{qj,X){yj_i - yj)] , 

0^ ^^ <^sNc , 
^ [QjA) = 

Since the t-channel factorised matrix elements are very fast to evaluate and the regularisation 
procedure does not add any complexity (because of the simple IR structure of the t-channel 
factorised matrix elements), the radiative corrections to all orders can be constructed as an 
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explicit phase space integral over each number of gluons emittecH: 

<^,Q^2j=Yi n(/ 



— , -- VJp ,=0 (2^)' J 2 J 

,=2 ^=l V-^P.x ^23) 



X X 



afA,q{Xa,Qa) X2fB,Qixb,Qb) (27r)'' dH ^Pk± C'2j({pi}). 



\k=l 



The cancellation of the poles in e ensures that the logarithmic dependence on A generated by 
the effective lower limit on the transverse momentum integrals cancels with the logarithmic A- 
dependence of the virtual + unresolved-real correction, which generates the exponential factors 
of Eq. ()19p . This is similar to the explicit construction of the solution to the BFKL evolution, 
where the very weak dependence of the solution on A at leading logarithmic accuracy was 
studied in Ref. [28p29]. and in Ref. [T9] at next-to- leading logarithmic accuracy. In Appendix [B] 
we investigate the stability under variations in A of a few of the cross-sections and distributions 
discussed throughout this paper. We find that the residual A-dependence is very weak — see 
Sec. [B] for further details. We then generally choose to use A = 0.5 GeV. Note that these 
findings are in good agreement with the conclusions from the studies of the A-dependence of 
the explicit solutions to the BFKL equations [19l[28l[29ll3T] , where the convergence of the phase 



space integration could be checked explicitly against an analytic solution. 

2 

The only remaining unregulated divergences of lA^jjljl related to the region of zero 
transverse momenta of the partons extremal in rapiditjQ. A similar situation was discussed 
in Ref. j21] . where simply a cut on the transverse momentum of the extremal partons was 
introduced, and the dependence of the cross section on this cut studied. We have refined the 
treatment for the current study. If there is no hard jet associated with the extremal partons, they 
could be viewed as not participating in the proper hard scattering of the event. In the parton 
shower picture, such emissions would be counted as (in this case) initial state radiation, and 
the divergence regulated by the Sudakov form factors. The treatment of these are beyond the 
scope of the current paper, and we will simply require that the extremal partons are associated 
with (i.e. a member of) a hard jet. With this requirement, the dependence on a lower cut- 
off of the momentum allowed for the extremal partons is weak. This is illustrated in Fig. [2] 
for a dijet-sample at the 7 TeV LHC, requiring justH two anti-kt-jets with absolute rapidities 
less than 4.5, and with transverse momenta above 30 GeV. We see that the contribution from 
transverse momenta much smaller than the jet scale is small. The requirement that the extremal 
partons be associated with a hard jet has to a large extent regulated the divergence for vanishing 
transverse momentum of the extremal partons (compare with Fig. 16 of Ref. [21j). In the results 
discussed in Sec. [5l we will choose a lower limit on the transverse momentum of the extremal 
partons, which is 5 GeV smaller than the minimum transverse momentum required on hard 
jets. Removing the very small contribution from smaller scales simply improves the phase space 
integration. Furthermore, the unregulated divergence at zero transverse momentum has to be 
explicitly removed. 



■^The lower limit on the transverse momentum in the phase space integrals is understood to be small, but 
non-zero, so Eq. (|2ip can still be evaluated numerically. 

^Actually, with the emission vertex of Eq. H10|l there is also a coUinear divergence for emissions close to the 
extremal partons from parts symmetrising pA ^ Pi and pB <^ Pn- We avoid this divergence by not averaging 
over the two contributions for emissions which are clustered into the same jets as the extremal partons. 

^Note that such a simple cut is problematic for NLO studies, because the truncation of the perturbative series 
introduces a large logarithmic dependence on any difference in the value of the cut applied on the two jets |42] . 
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Figure 2: The dependence of the dijet cross section on the transverse momentum of the parton 
of lowest or highest rapidity. The step at 30GeV is caused by the requirement of minimum 
30GeV transverse momentum by the hard jets. 

The construction of an exphcit integration over emissions to all orders relies on an efficient 
phase-space generator |2Uil21j . which should sample final states with the number of particles 
varying by more than an order of magnitude. The problem is significantly different to that 
of a so-called general purpose Monte Carlo (Pythia [IS], Herwig [H], SHERPA [45j), since in 
these approaches, the approximation to the virtual corrections is defined such that the emission 
of particles is unitary, i.e. does not change the total cross section, which allows for a simple 
probabilistic interpretation. In HEJ, an approximation to the virtual corrections is calculated, 
and introduces a suppression of the regularised matrix element for any final state with a finite 
number of partons, as the rapidity length of the event is increased. This is countered by the 
(positive) contribution from the emission of additional gluons, and introduces a correlation 
between the number of final state partons and the typical rapidity length of an event. It is 
absolutely crucial to incorporate this probabilistic correlation in the phase space generator in 
order to obtain satisfactory numerical stability in a finite amount of time. Such a phase space 
integrator can be efficiently implemented by following the ideas of Ref. [40j . The fully exclusive 
formulation in a flexible Monte Carlo facilitates the study of any observable. 

3 Matching 

The previous sections have set up the all-order approximations to jet production of High Energy 
Jets, and discussed the implementation as a flexible Monte Carlo, integrating explicitly over 
n-particle phase space. The resummation procedure generates only certain partonic phase space 
configurations (FifL-configurations, see Section [2.2.ip . The dijet production process is calculated 
within this approximation as (for notational brevity, we have omitted the label indicating the 
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use of the regularised amplitudes) 



/^J2 n=2 i=l ^•^PiX=0 (^^-^ J J ^ 

afAJi{Xa,Qa) a;2/B,/2 (Xfe, Qfe) (27r)'' (J^^^Pfc±^ C'2j({Pi}), 



(24) 



X Xo 



where the first sum is over the flavours /i,/2 of incoming partons. The distribution of any 
observable can be obtained by simply binning the cross section in Eq. (|24p in the appropriate 
variable formed from the explicit momenta. Obviously, multi-jet rates can also be calculated by 
multiplying by further multi-jet observables O^j, O^j, ... in Eq. (p^ . 

In section [3TT] we will discuss how the amplitudes for the FKL-states included in Eq. ()24p can 
be corrected to full tree-level accuracy, limited only by the availability of full tree-level matrix 
elements. In Section 13.21 we will discuss the inclusion of all remaining partonic configurations 
(in practice for up to 4 jets). 



3.1 Matching for FKL Configurations 

Firstly, we want to match the description of the FKL n-jet configurations to the full tree-level 
matrix elements and thus improve upon the approximations inherent to the resummation. This 
can be straightforwardly done because of the flexibility inherent in Eq. (|24p . Let 0^j{{pi}) 
denote the measurement function for exclusive n-jet production acting on the partonic phase 
space. This function will return one if the chosen jet-algorithm finds exactly n hard jets in the 
m-partonic phase space point, and returns zero otherwise. Furthermore, it will give access to the 
momenta of the n jets, {pjii{Pi})}- We note that 02j{{pi}) = Y^'^=2^nji{Pi})- principle, 
we would then want to simply multiply each exclusive jet measure function with 

\Mhh^h^-9fmpj^{{p,})})f ^^^^ 
|A^t,/i/2^/i<;-9/2({p^^({p,^})})f 

where the numerator is simply the (spin and colour summed and averaged) square of the full n- 
jet tree-level matrix element, and the denominator is the i/ii^J-approximation to this tree-level. 
This would ensure tree-level accuracy of the n-jet rates, while simultaneously weighing the n-jet 
samples with the virtual corrections from HEJ. 

However, a few modifications to this naive approach are necessary. Firstly, the jet momenta 
may not be of zero invariant mass. Secondly, the transverse momenta of the jets generally 
will not sum to zero, since some of the partons generated in the event may not be included in 
the hard jets. We therefore have to construct a new set of n jet-momenta to be used in the 
matching. We start by making each jet momentum equal to the sum of the parton momenta of 
each jet (each jet contains mostly just one hard parton after the ilTE'J-resummation) . We then 
redistribute the transverse momenta of any partons not belonging to a jet among the hard jets, 
and remove these softer partons from the list of particles (and momenta) used in the matching. 
We choose to distribute the momenta in proportion to the transverse momenta of the resolved 
jets. If the sum of the momenta of the non-jet partons is q and the scalar sum of the transverse 
momenta of the jets is P±, the new set of hard momenta p^^ is given by 

pT=PJi+'1*^-^- (26) 



15 



The energy component of each jet is then finally reset to put it on-shell, and the momenta of 
the incoming partons are defined by energy /momentum conservation. 

This reshuffling of momenta is illustrated for a sample event in figure [3l which has eleven 
partons in the final state, in a momentum configuration leading to four hard jets with transverse 
momentum above 30 GeV, found with the anti-kt jet algorithm, as implemented in Fast Jet [39] . 
The red circles show the positions in rapidity-phi space of the partons; the radii of the circles are 
proportional to the transverse energy of each parton and jet (and do not, therefore, represent 
the area of each jet). The green circles indicate the jets of the original event. As expected, 
they coincide with the hardest quarks/gluons. The blue circles indicate the reshuffled momenta 
used in the matching. Note, this procedure does not change the kinematics of the actual event; 
only the reweighing of the event to full tree-level accuracy is performed with matrix elements 
evaluated for the slightly modified momenta. If the threshold on the transverse momenta of 
jets was set very low, and the jets were finely resolved (small i?-parameter), then no reshuffling 
of momenta would be necessary. However, the full matrix elements can only be evaluated for 
states of relatively low multiplicity (with MadGraph [46] , we limit ourselves to matching of up 
to four jets). So with a low jet matching scale, the available fixed order matrix elements for 
matching would cover only a small part of the total cross section. A similar issue occurs for the 
CKKW-L [6ll3 or MLM jS] style matching of parton shower algorithms. 

We then reweigh each event generated with the following multiplicative matching factor, 
evaluated with the on-shell hard momenta as found by the described procedure: 



Wn-iet = ^ - (27) 



In this notation, we have suppressed the flavour and momentum-dependence of Wn, but it is 
obviously calculated on an event-by-event basis. The FKL-matched cross section is then found 
as 



" V4,=A (2-)3 i l) 



E En 

/i,/2 n=2 i=\ 



resum, match \ ^ \ ^ 1 T / / CI P},J_ 



X 

m 



X X 



^Ol,j{{pi\) w^--^^^ (28) 

1 

afA,h(.Xa,Qa) X2/iJ,/2 (Xfe, Qfo) (27r)^ '^Pi± C>2j(.{Pi})- 



The impact of this matching procedure can be seen in Fig. [H which displays the differential 
dijet cross section wrt. the rapidity difference Ayfi, between the most forward/backward hard 
jet, within the following set of cuts: 

Pj^ > 60GeV \yj\ < 4.5 anti - kt, i? = 0.6. (29) 

The matching scale is set equal to the general jet scale of 60 GeV. The red (dot-dashed) curve 
is the result of the pure resummation; the blue (dashed) curve is obtained after matching of 
the states arising in the resummation up to four hard jets. The correction is small throughout, 
being slightly more significant at low rapidity spans. 
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Figure 3: This plot shows for an example event the momenta of the partons (red), the resultant 
jets from Fast Jet (green) and the reshuffled momenta described in the text (blue). The radii of 
the circles are proportional to the transverse momentum of the particle or jet described. 




Figure 4: This plot shows the impact of matching in jet production, as a function of rapidity 
span. The initial HEJ approximation is shown (red, dot-dashed) together with matching to FKL 
configurations (blue, dashed) and matching to both FKL and non-FKL configurations (green, 
solid). 
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3.2 Matching for Non-FKL Configurations 



The processes and partonic configurations whicli do not arise in the resummation are included 
straightforwardly by adding these to the dijet rate found by the (matched) Eq. (p^ . For example, 
we can add the remaining contribution to the exclusive dijet rate as 



2 



a. 



non-FKL _ \ " \ ^ 

2i 11 W (2vr)3 J 2 

/l,/2//l,//2 *=1 '^•'P^-L-P-Lmm V j 



- $B /^) l-^"-^:f ,30) 



a; Qa) X2fBJ,{Xb,Qb) (2^)" ^^\JlP^^j (^2j{{Pi}), 

where p_Linin is the minimum transverse momentum required for hard jets. The function B({/j}, {pi}) 
returns one if the parton and momentum configuration is of non-FKL status. If only rapidity 
ordered sets of momenta pi are generated, then one needs to also sum (or Monte Carlo sample) 
over all possible assignments between momenta and the particles in the process. The generali- 
sation to the three and four jet states is straightforward, and the final result for the dijet rate 
is 

resum, match , \ ^ non~FKL /"oiN 

= + 2^ a^j . (31) 

n 

Each component is implemented by explicit Monte Carlo sampling over phase space and an 
evaluation of matrix elements. Therefore, any observable can be constructed and studied, also 
after matching has been included in the formalism. 

The impact of the non-FKL states is indicated on Fig. [H where the green (solid) line is 
obtained from the sum of all terms in Eq. ()3ip . This correction is again more significant for 
small rapidity spans, as we expect. 



4 Logarithmic Corrections to the Scale Choice 

The discussions so far have made no assumptions on the scale choice made for the evaluation 
of as or the pdfs. In this section we will compare the results arising for a fixed scale choice 
(of e.g. the minimum jet transverse momentum), and a scale choice made event by event equal 
to the maximum jet transverse momentum of the event. Finally, we will include pieces of the 
next-to-leading logarithmic corrections to the BFKL kernel, which will stabilise the dependence 
on the scale choice. This will then form the basis of the standard scale choice for the results 
presented in Section [S] 

The connection between the formalism of High Energy Jets and that of BFKL |14ftl6p33ll47j is 
that in the limit of large invariant mass between all partons (conditions relaxed for neighbouring 
pairs of particles at NLL [48j), then the amplitudes underlying the BFKL formalism coincide 
with those of HE J (and with those of full QCD). The NLL corrections to the BFKL kernel 
have two origins: the one-loop corrections to one-gluon emission, and the contribution from 
two-gluon and quark-anti-quark-emission in quasi-multi-Regge-kinematics (i.e. not necessarily 
a large invariant mass between the pair of particles). The net result of the corrections is a 
sum of an expression with the same functional form as the LL kernel, multiplied by a running 
coupling logarithm, and a term of a more complicated kinematic structure [18pi9j . The relevant 
discussion of the regularisation of the NLL corrections to the BFKL kernel was presented in 
Ref. |18lll9t BT]. We repeat it here, with a notation tailored to the present application. 
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The NLL BFKL kernel is expressed in terms of a transverse momentum, which is the trans- 
verse momentum of the emitted gluon, or the sum of transverse momenta of the emitted pair 
of gluons or quark-anti-quark pair. All other kinematic dependence is integrated over before 
arriving at the BFKL kernel. In D = 4 + 2e the BFKL amplitudes obey the following relation 
at NLL accuracy (compare with Eq. (|17p ) 



A/JBFKL 

Pa Pb^Pl 



Pi-1 Pi Pi+1 ■■■ Pr, 



A/fBFKL 

'Pa Pb - 



Pl ••• Pi_l Pi + l 



(32) 



with lC{pi) = J^liPi) + where K'^{pi) is irrelevant for the current discussion, and 
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(33) 



with /3o = ^-^c ~ I'V- '^^^ NLL-corrections to the trajectory give 



a(g^ 



(/3o/iVc 



/vr2 67\ /202 



vr 



18 



C(3) e 



(34) 



3Afc 



28 



vr 
"3 



By applying the same regularisation procedure as discussed in Sec. 12.41 we find that for both the 
real emission (evaluated above p? > A^) and for the trajectory the term found at LL accuracy 
is multiplied by a running coupling logarithm. For the real emission this is: 



(4vr) 4vr 



(35) 



For the regularised trajectory we find 

as Ca 



vr 



ln(^ 



2 4vr q^X^ 



(36) 



These results are in complete agreement with what was found in Ref. [18>il9j . The logarithm 
of the trajectory may seem a little odd (being dependent on A), but it reproduces the NLL 
BFKL results when expanded in /3o- Besides, the study of the pure NLL BFKL correction in 
Ref. [18 1 119 ^ ^9] , show that the organisation of the cancellation of soft divergence is completely 
stable for the values explored for A. 

We finish off this section with a simple study illustrating the impact of various scale choices on 
the average number of hard jets versus the rapidity difference between the most forward/backward 
jet within the cuts of Eq. ()29p . We apply three different choices: 1) a fixed scale choice of 60GeV, 
2) a common scale choice, chosen event by event, of the largest transverse momentum of any 
jet, and 3) the latter, including the logarithmic corrections discussed above. 

This observable is just one of many with a strong correlation with the number of hard jets - 
in order to describe the region of phase space of large Ayfh, it is clearly imperative to describe 
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correctly the emissions of many hard jets. Other such examples are studied in the next section. 
The obvious expectation is that for a larger value of as (smaller scale), one would see more hard 
jets than for a smaller value of Os (larger scale). Indeed, this is found in Fig. [5] for rapidity 
differences less than roughly 5. Furthermore, we see that including the logarithmic corrections 
outlined above leads to a prediction in-between that of the fixed, low scale choice of 60GeV, and 
the choice of the hardest jet scale. This is of course entirely as expected. For larger rapidities, 
phase space constraints become increasingly important, and the scale at which the pdfs are 
evaluated will influence the details. Eventually, as Ayj^ increases further, the average number 
of hard jets decreases as the phase space for additional radiation is reduced when the energy of 
the forward/backward jets gets close to the total available hadronic energy. 

5 Results 

In this section we present results for dijet- and trijet-studies within these cuts (identical to those 
of Eq. ([29]) in Section [3T]) : 

Pj^ > 60GeV \yj\ < 4.5 anti -kt,R = 0.6. (37) 

We choose as the standard scale choice fir = maxj(p^j), and include the running coupling 
logarithms from NLL accuracy, as discussed in Section HI As shown explicitly in Appendix[Cl the 
all-order framework of High Energy Jets is free from the instability seen in the NLO-calculation 
of dijet production [32] when the transverse momentum cut on the two jets is equal. This 
problem simply arises from the fact that in a three-particle system, the cut on the transverse 
momentum of two particles automatically changes the phase space explored by the transverse 
momentum of the third particle. The infra-red region of the real emission corrections to the dijet 
system is explored in the limit where the two transverse momenta of the hard jets are equal. 
An off-set Ap± in the cut of the two hardest jets modifies the soft phase space for additional 
real emission, and can therefore introduce a logarithmic dependence on Ap±. However, this 
dependence seems specific to cross sections terminated at NLO, and is washed away in several 
other all-order frameworks, e.g. POWHEG |10p50p51| (NLO matched to a parton shower) and the 
BFKL generator studied in Ref. [30J. Since the problem is related to a fixed-order perturbative 
calculation rather than any observation or our description, we will proceed with an equal cut on 
the transverse momentum of all jets. 

We will apply the anti-kt jet-clustering algorithm as defined and implemented in Ref. [32] 
with R = 0.6. We define Ayfb as the rapidity difference between the most forward and most 
backward hard jet. The average number of jets in the events is an obvious indication of the 
importance of the hard, higher order corrections that are resummed in High Energy Jets. In 
Section 15.1.11 we study simple characteristics of the inclusive sample generated with HE J; we 
then move on to discuss distributions in p±, Hj- (scalar sum of transverse momenta) and Sij 
(invariant mass between hardest jets), where the corrections have a particularly large impact. 
Other all-order approaches like e.g. Cascade [52^53] calculate higher order corrections in the kt- 
factorisation scheme through the evolution of off-shell pdfs convoluted with a 2 — )• 2 (off-shell) 
hard scattering matrix element. It would be interesting to compare the predictions for these 
observables also from such a framework. 



20 



(0 
0) 

3 

o 



' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

Inclusive Dijet Production 
Vs=7 TeV, anti-kt(R=0.6) 
p.>60 GeV, |y.|<4.5 



--■H_r=60GeV 

Ji_f=max(p^|) 
H_r=max(p^|)*log 



2.8 



2.6 



2.4 



2.2 




Figure 5: The average number of hard jets {p± > 60 GeV) in inclusive dijet production as a 
function of the difference in rapidity Ayfb between the most forward/backward hard jet. 




Figure 6: The differential dijet-rate with respect to the rapidity difference between the most 
forward and most backward hard jet (left) and the transverse momentum of the hardest and 
second-hardest jet of the event (right). 
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5.1 DiJet Studies 

5.1.1 Rapidity and Transverse Momentum Distributions 

In Fig. [6] we have plotted the differential cross section with respect to both Ayj^, (left) and 
the transverse momentum of the hardest and second-hardest jet in the event (right). The dijet 
rate is peaked at zero rapidity difference, and the radiative corrections have clearly induced a 
difference in the transverse momentum distribution of the hardest and second-hardest jet (which 
is obviously identical at leading order). The transverse momentum spectrum is compared to that 
arising in a LO calculation (using the MSTW2008LO pdf set, and setting the renormalisation 
and factorisation scale equal to p±j)- The LO spectrum is significantly softer then that of the 
hardest jet arising in HE J. 

5.1.2 ^Vfbi Ht, Sj-^^j2 and the Average Number of Jets 

In Fig. [5] we plotted the average number of hard jets (transverse momentum larger than 60GeV) 
according to the anti-kt jet algorithm with R = 0.6 in the inclusive dijet sample, as a function 
of the rapidity span Ayji, between the most forward and most backward hard jet. As expected, 
there is a strong correlation between Ayf^ and the average number of hard jets. The average 
number of hard jets rises monotonously until Ayfb ~ 7, simply because the partonic phase space 
increases. However, as the rapidity span is increased further, the parton density functions fall 
off so steeply as x — )• 1 that the production of additional hard jets beyond the required dijet 
system is effectively vetoed. 

One observable which is often used in the search for signals of new physics at hadron colliders 
is the scalar sum of transverse energy (or momentum) in the hard event. For the jet studies, we 
define it as 



where the sum runs over the jets found with a given jet-algorithm, with a transverse momentum 
bigger than some hard cut-off. In Fig. [7] (top right) we plot the differential cross section wrt. Ht 
as obtained both at leading order QCD, and within HEJ. The distribution is clearly more 
pronounced at large Ht when the higher order corrections from HEJ are included. This is made 
very clear on the plots of the normalised i?r-distribution at the top right of Fig. [71 The bottom 
plot in Fig. [7] is of the average number of jets in the events as a function of Ht- We see that the 
average number of jets starts at 2, and very quickly rises above 3 (already at roughly GOOGeV). 
A priori, one might have expected the large- -ffj tail to be dominated by two hard jets. Fig. [7] 
clearly demonstrates this is not the case. Furthermore, the very high average number of jets in 
the large- f^T tail of the dijet distribution suggests that a veto on further hard jets beyond two 
would be very efficient in suppressing the QCD contribution to large- -ffr dijet events. The rise in 
the number of hard jets is a direct consequence of the i-channel colour exchange, and therefore 
may be different between the QCD process and any process originating from new physics. 

In Fig. [8] we plot the same three quantities for Sij, the square of the invariant mass between 
the two hardest jets of the event. From the top-right plot of the normalised distribution we see 
that the corrections implemented in HEJ lead to a relative enhancement at small Sij, and a 
suppression at large Sij compared to the LO result. This behaviour can be explained by the fact 
that the radiative corrections implemented in HEJ will fill the rapidity span between the two LO 
jets, as seen in e.g. Fig. O The hardest jets are likely to be the ones that are radiated centrally in 
rapidity (a similar effect was seen in Ref. [21j). For any value of the rapidity span Ayjf, between 
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Figure 7: The differential cross section witli respect to Ht (top left), and the normalised spec- 
trum (top right). The radiative corrections implemented in High Energy Jets enhance the 
high-Z/j'-tail significantly. The bottom plot is of the average number of hard jets (transverse 
momentum above 60GeV) in the events as a function of Ht- Hard radiative corrections, as those 
included in HE J, are clearly important in the description of events with large Ht- 
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Figure 8: The differential cross section witli respect to the square of the invariant mass Sij 
between the two hardest jets (top left), and the normalised spectrum (top right). The bottom 
plot is of the average number of hard jets (transverse momentum above 60GeV) in the events 
as a function of Sj,. Please see text for discussion. 
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the most forward and most backward hard jet, the higher order corrections implemented in HE J 
can produce a central jet, which has a slightly harder transverse momentum spectrum than the 
extremal ones, and is therefore more likely to be the one(s) used in the definition of Sij. For 
all values of Ayj^j the iJE'J-corrections will lead to a smaller value Ay^ of the rapidity span 
between the two hardest jets in the event. This in turn induces a smaller value of Sij than would 
be the case in the description of LO exclusive dijets. 

The average number of jets versus Sij is shown at the bottom of Fig. [8l At small values 
of Sij it is peaked at almost 3.8, then falls off abruptly to 2 already at Sij = 2 • lO^GeV^, and 
then rises to a plateau at 2.25. This behaviour is a sum of two effects. The events at small 
Sij are dominated by the cases where two additional (central) jets have been radiated by the 
HE J -mechanism; these two jets are often the hardest (in pt), and therefore define Sij (which 
can be small since the jets can be aligned in p± (i.e. they do not have to be back-to-back in 
azimuth as dijets at LO) and close in rapidity). The strong peak at small Sij is therefore an (at 
least) a^-correction to the tree-level dijets. There is another effect, giving rise to a distribution 
increasing with Sij, starting at 2 for Sij = and then reaching the plateau. This is just the 
standard C'(as)-effect of one hard radiation. The value 2.25 is not too far from the value for the 
average number of jets in the fully inclusive dijet sample generated with HEJ, so the plateau 
is just an indication of only a small correlation between Sij and the average number of jets. 
However, for the C'(as)-correction of additional jets to arise and the average number of jets to 
rise from the LO value of two, a rapidity difference between the two jets is required (c.f. Fig. [5]), 
and this naturally leads to an increase in Sij. This why the average number of jets is close to 2 
only for small Sij. 

5.2 Trijet studies 

In this section we will briefly present the distributions discussed previously, but this time for 
events with at least three hard jets. Many of the features identified in inclusive dijet production, 
like e.g. the strong correlation between the average number of jets and the rapidity difference 
between the most forward and most backward hard jet are found also for 3-jet production. 

On the top of Fig. [9] we compare the results for the i^T-distribution in 3-jet production 
as obtained in leading order and in High Energy Jets. Top left is the distribution in absolute 
numbers, on the top right the distribution is divided by the total 3-jet cross section. Just as in 
the dijet case, the higher order corrections implemented in HEJ hardens the spectrum in Ht- 
The average number of hard jets vs. Ht is shown on the lower left of Fig. [9l The distribution 
rises to 3.6 at Ht = 800 GeV and then drops off slightly for increasing Ht- We note that the in 
the trijet case, the average number of jets rises 0.6 units above the minimum required, whereas 
in the dijet case it rises a full 1.2 units. Both cases represent of course large corrections to the 
simplistic tree-level point of view. The lower right plot on Fig. [9] is of the average number of jets 
vs. the rapidity span between the most forward/backward hard jets. It rises from 3 to roughly 
3.8 at rapidity differences between 7 and 8, before dropping back down towards 3, again because 
the increase in x necessary for additional radiation leads to a pdf suppression. An increase in 
the centre-of-mass energy of the proton-proton collision will obviously lead to a further increase 
in the number of hard jets produced. 

In Fig. [To] we study the distributions in the invariant mass between the two hardest jets in 
the event. Similarly to the dijet-case, we find that the results from HEJ are suppressed at large 
Sij compared to the 3-jet LO estimate. The bottom plot on Fig. [10] is of the average number of 
jets in inclusive 3-jet events as a function of the invariant mass between the two hardest jets. 
The distribution is very strongly peaked at small Sij for exactly the same reason as the dijet 
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Figure 9: The differential 3-jet cross section with respect to Ht (top left), and the normalised 
spectrum (top right). The radiative corrections implemented in High Energy Jets enhance the 
high-Z/r-tail significantly. The bottom plot is of the average number of hard jets (transverse 
momentum above 60GeV) in the events as a function of Ht- Hard radiative corrections, as those 
included in HE J, are clearly important in the description of events with large Ht- 
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Figure 10: The differential 3-jet cross section with respect to the square of the invariant mass Sij 
between the two hardest jets (top left), and the normalised spectrum (top right). The bottom 
plot is of the average number of hard jets (transverse momentum above 60GeV) in the events 



as a function of Sij. Please see text for discussion. 
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case: the correction from additional jet production allows for two central (and thus generally 
slightly harder in transverse momentum) jets, which can form a system of small invariant mass. 
In fact, the average number of jets in inclusive 3-jet production at small Sij is almost identical 
to the average number of jets in inclusive dijet production at small Sij. 

5.3 Gaps Between Di-Jets 

It is possible to construct several observables which are sensitive to additional radiation from 
the dijet-system, and thus can be used as a direct test of the description arrived at using various 
descriptions like e.g. fixed order, shower. Cascade |52ll53j and HEJ. There is a small challenge 
in defining quantities which are stable within each perturbative framework. We have already 
discussed that the NLO calculation for dijet production is unstable in a setup of equal transverse 
momentum cuts on the two jets [30lll2], but that it can be stabilised by requiring e.g. cuts of 
65 GeV and 60 GeV on the hardest and next-to-hardest jet. Some dijet observables will then 
be calculable in several frameworks. A well-studied example of an "inclusive" dijet observable 
is the average of cos(7r — vs. Ayfh. At LO, the two jets are obviously back-to-back, and 
cos(7r — (pjj) = 1. The benefit of this observable is that it is completely inclusive in the radiation 
between the two jets (this inclusiveness allows for studies also within (semi-) analytic approaches 
of BFKL |54^l55j) - any emission will cause a decor relation, whether or not it is identified as 
a separate jet. An experimental study of this quantity could serve as a strong test of the 
description of higher order corrections. 

Instead of studying the effect of additional radiation through its impact on the jets extremal 
in rapidity, one can study the radiation in-between the dijets directly. The Atlas collaboration 
have published a note on such a study in early data from the LHC |56j . They present data for 
the so-called "gap fraction" , defined as the fraction of dijet events with no additional hard jets 
between the two (in rapidity). We have already seen (e.g. Fig. [5]) that the average number of 
hard jets in dijet events increases with the rapidity difference between the forward/backward 
jet, and this should obviously be reflected in the "gap fraction". These early studies also serve 
to guide jet veto studies for Higgs boson production in association with dijets OEZ]. For these 
studies, it is of interest to use a small transverse scale for the vetoing of further jets. The Atlas 
study defined jets using the anti-kt algorithm, with R = 0.6, and a transverse scale of 30 GeV. 
In order to ensure a sufficiently small dijet-rate (and thus an acceptable scaling factor for the 
trigger), a harder scale was required. We will here concentrate on the part of the study, where 
the average transverse momentum of the two jets extremal in rapidity was required to be above 
60 GeV. The cuts used are then: 

Pj^ > 30GeV p± > 60GeV \yj\ < 4.5 anti -kt,R = 0.6. (39) 

where pj_ is the average transverse momentum of the most forward/backward jets. 

Figure [11] presents the prediction for both the average number of hard jets (with a transverse 
momentum larger than 30 GeV) and the gap fraction obtained using HEJ, within the cuts in 
eq. pop . We have also indicated the variation in both quantities between a scale choice of 
30 GeV, max(p^j) and of max{p±j) including the logarithmic corrections discussed in SecHl all 
using the pdfs included in MSTW2008 [58j. For the last, "central" scale choice we also present 
the results obtained by using NNPDF2.0 [59j, including the full envelope of the 100 uncertainty 
pdfs. The uncertainty induced by the pdfs on these quantities is completely negligible (they 
begin to play a role at Ayfj, > 8). The uncertainty estimate induced by a variation in the 
renormalisation and factorisation scale between 30 GeV (the minimum transverse scale for jets) 
and max(p_Lj) is increasing for increasing rapidity spans, and amounts to a variation between 
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Figure 11: The average number of jets and gap fraction vs. the rapidity difference Ayj^ between 
the most forward and most backward jets. The upper and lower sohd Unes are for scale choices 
of 30 GeV and m.ax(p±j) respectively. The central solid lines are for a scale choice of m.ax{pj_j) 
plus the logarithms of Sec. H] for MSTW2008 [58j uncertainty pdf sets. The difference between 
the sets is barely observable. Also shown in dotted and dashed lines are the results with the 
uncertainty sets from NNPDF2.0 [59{ . 

3.6 and 4.0 in the average number of jets (with a transverse momentum larger than 30 GeV) 
for a rapidity span of 7. These results are marked by the outer solid lines (on both plots). The 
central solid line is obtained by choosing the renormalisation scale m.ax(p±j), but including the 
logarithms as discussed in Sec. [H The results obtained by choosing the renormalisation scale as 
30 GeV and including the logarithms is almost identical. 

The pdf and scale uncertainty of the predictions for the average number of jets and the 
gap fraction are sufficiently small that the the ideas and calculations presented here can be 
meaningfully confronted with data, once it has become available. 

6 Summary and Conclusions 

We have discussed the implementation of the framework of High Energy Jets |lipi2 ] in a flexible 
Monte Carlo; the new components discussed in the present paper include 1) the organisation of 
the all-order cancellation of IR divergences between real and virtual corrections (Sec. I2.4p . 2) 
matching to high multiplicity tree-level matrix elements (Sec. [3]), and 3) the inclusion of higher 
order logarithmic terms to stabilise the scale dependence (Sec. H]). 

In Sec. Owe studied the impact of the perturbative corrections included in High Energy Jets 
on a number of dijet and trijet distributions. We find that compared to LO, the distribution is 
harder in both the transverse momentum and in Ht (the scalar sum of transverse jet momenta), 
while the invariant-mass distribution between the two hardest jets is softened. Similar results 
hold for trijet-observables. Therefore, the understanding of the radiative corrections could lead 
to better methods for suppressing the Standard Model contribution in new-physics searches. 

The effect of hard, perturbative corrections is cleanly displayed by the average number of 
hard jets versus the observable in question. This is particularly true for the rapidity span 
between the most forward/backward hard jets, which is a direct measure of the phase space 
available for hard radiation. 

Finally, in Sec. 15.31 presented the prediction obtained from High Energy Jets of an 
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observable sensitive to inter-jet radiation, which is currently under study by Atlas [56]. We find 
that the theoretical uncertainty on the quantities studied is dominated by the scale choice, while 
the variation induced by pdf uncertainties is completely negligible. 

The generator High Energy Jets is available at http://cerii.ch/hej. 
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A The Building Blocks for the Regularised High Energy Jets- 
Cross Sections 

We define here in one place the necessary building blocks to construct an amplitude in the HEJ 
framework: 



|-^HEj(ta})r = 4 (jyj _ \\Shf2^hf: 



rcg - . . „ . I|2 



n— 2 / / 1 A \ \ f4n'> 

■ n [s'^Ca (^-^V^^{q^,q^+i)V^{qi,q^+i) - ^ 9 (pf < X^)jj 

n-l 

• Jl exp [uj^iqj, \){yj-i - Vj)] . 
j=i 

Our momentum convention will be that pA and pB represent the momenta of the forward and 
backward moving initial partons respectively. The outgoing momenta of all quarks and gluons 
are then numbered in decreasing rapidity so pi is the most forward etc. We then define qi to 
be the momenta which correspond to the t-channel momenta in the effective t-channel exchange 
picture, that is 



qi=PA-Pi, qi = qi-i-pi 2<i<n-l. (41) 

nt pieces, 

quarks this is 



The current pieces, ll'S'/j/j^i./j/jH^, indicate the square of pure current-current scattering. For 



II c l|2 _ I — — 1*^ _i_ I — •+ 1"^ _i_ I •+ — 1*^ _i_ I •+ •+ 1^ //lo^ 

W^qQ^qQW —\Jal'hri\ ^yal'hri] ^yal'hri] '^Wal'hn\ ' V^-^J 

Anti-quarks are treated in the same way with — )• = v~^^^v^. For gluons it is more 
complicated as there is an overall factor for helicity conserving channels (see eq. ([9])) and also 
flipped-helicity contributions weighted with different factors |12j : 

WSqg^q^? = (2 (1 - ) + 

(\- -|2 , I- -+12 , !•+ -|2 , |-+ (43) 
\\3al'hn\ ^yal'hrA ^\3al'hn\ + Mai " Jfon I j 
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For this example of a backward- moving parton z = p^/p^ ■ The result for a forward-moving 
gluon is the same with Pa Pb, Pi Pn and z = pf /p^. The helicity- flipped contributions are 
given by 



\S, 



|2 1 



- 1 
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(45) 



Gluon-gluon scattering is the natural generalisation of what has gone before. There are now two 
relevant ratios zi = pf /p^ and Z2 = P2 /Pb • Then we define 
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(46) 



As the necessary factor has already been included explicitly in eq. (j46p . the helicity-specific terms 
I 1 2 

like \Sg-g-^g-g+\ arc equal to the corresponding quark-gluon ones if the helicity-conserving 
gluon was replaced by a quark, e.g. 
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9 9 -^9 g^ 
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(47) 



We do not include the contributions where neither gluon conserves helicity (although this could 
be generalised) because the contribution is so insignificant. 
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Returning to the remaining pieces of eq. (j40p . the factors Kf^ are straight-forward and 
inspired by the exact high-energy hmit: Kg = Cp for a quark of any flavour and Kg = Ca for 
gluons. 

The emission vertices were given in eq. (|lUp : 

V^{qi,qi+i,PA,PB,Pi,Pn) = - {qi + qi+iY 

2 \Pi+l-PA PA-PB PA-Pn J (4»j 

Pb ( , • PA , • Pi \ 

' H \ ]-PB^Pn- 



2 \Pi+i -pB PB- PA PB ■ Pi 

The combination of the virtual corrections in Section [XHl and the regularisation in Section [X 
give the following factor in the exponential in eq. ()40p : 

-^.„A)= -^log| (49) 

where the bold indicates that it is the sum of the square of the transverse components which 
are included in the log. 



B Variations of the regularisation parameter A 

In this appendix, we show a few results to demonstrate that our conclusions are not sensitive to 
the chosen value of the regularisation parameter A. This is the scale above which radiation is 
considered to be a real emission. The regularisation procedure is described in full in Section [2.41 
We again use the cuts defined in eq. ([29|) throughout. 

Table [T] shows the exclusive n-jet rates for n = 2, 3, 4, for values of A from 0.2-2 GeV. We 
can see that the changes in A do not have a large effect, particularly for the 2j and 3j rates. 



A (GeV) 


a(2j) (/ib) 


a(3j) (/ib) 


a(4j) (/.b) 


0.2 


1.58 


5.90 E-2 


9.6±0.1 E-3 


0.5 


1.58 


5.93 E-2 


lO.libO.l E-3 


1.0 


1.59 


5.95 E-2 


9.7±0.2 E-3 


2.0 


1.61 


5.99 E-2 


9.2±0.2 E-3 



Table 1: Exclusive n-jet cross sections for different values of the regularisation parameter A. 
The errors shown are statistical - they are not shown for the 2j and 3j rates because they are 
smaller than the last quoted digit. 

Figure [12] shows the distribution of the rapidity difference between most forward and most 
backward jet, Ayfj,, for different values of A both for the inclusive 2-jet sample and the exclusive 
3-jet sample. The differences are very small. We use A = 0.5 GeV as the default. 

In the HEJ framework, the number of quarks and gluons is treated as a variable and contri- 
butions are summed over n from 2 to oo, see eq. ([23|) . In practice, there is an upper cut-off on 
the value of n. This has a very weak dependence on A as it stands to reason that the lower the 
cut-off on resolved emissions, the more of them you need to consider to get the same results. 
For A = 0.5 GeV, we use rimax = 22, and find no observable difference in physical results by 
varying around this value. 
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Figure 12: This plot shows the ^yfb distribution, for different values of the regularisation 
parameter A for the inclusive 2-jet sample and the exclusive 3-jet sample (times 40). The 
differences are small and we choose to use a default value of 0.5 GeV. 

C Stability of Equal Cut in the Transverse Momentum of Dijets 

As discussed in Section 15. H NLO calculations suffer from an instability when the transverse 
momentum cuts on the two jets are symmetric f42j . The same effect however is not seen in 
other calculations which are not terminated at NLO [301151]. In line with this, HE J does not 
see an instability here either. We illustrate this by implementing an offset-parameter, Ap_|_, and 
requiring a difference in the transverse momenta of the two jets of at least this value. Figure [131 
shows the dijet cross section as a function of Ap^. It continues to increase linearly as Ap_|_ 
decreases to zero and shows no sign of the turnover observed at NLO. 
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